The influence of gas-kinetic evolution on plasma reactions 
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Plasmas in which there is a threshold for a dominant reaction to take place (such as recombination 
or attachment) will have particle distributions that evolve as the reaction progresses. The form of the 
Boltzmann collision term in such a context will cause the distribution to drift from its initial form, 
and so cause for example temperature fluctuations in the plasma if the distribution is originally 
Maxwellian. This behaviour will impact on the relevant reaction rates in a feedback loop that 
is missing from simple chemical kinetic descriptions since the plasma cannot be considered to be 
isothermal, as is the case in the latter approach. In this article we present a simple kinetic model that 
captures these essential features, showing how cumulative differences in the instantaneous species 
levels can arise over the purely chemical kinetic description, with implications for process yields and 
efflciencies. 

PACS numbers: 52.20.Fs, 52.20.Hv, 52.25.Dg, 52.25.Gj, 52.25.Ya, 52.27.Gr, 82.33.Xj 



INTRODUCTION 

In many applications and contexts plasmas that evolve 
to a critical non-equilibrium stage are vitally important, 
particularly where such plasmas develop to a stage where 
either an instability is triggered (such as the onset of a set 
of reactions that alter the plasma characteristics) or a sig- 
nificant phase change is induced (the plasma evolves to a 
substantially different equilibrium condition). Examples 
of such behaviour are when plasmas become significantly 
electronegative [for example[l|, 0, H or quench completely 
to form neutral gas [for example |j, d, 0, 111 ■ The pur- 
pose of this article is to explore the implications in the 
gas-kinetic framework of plasmas in which the distribu- 
tion function of the electrons is altered by unbalanced 
plasma chemistry; reactions with strong electron energy 
dependence such as electron recombination with positive 
ions to form neutrals, or electron attachment to neutrals 
to form negative ions, can lead to significant changes in 
the electron energy distribution function (eedf) which in 
turn affect the progress of subsequent reactions. 

Very often the modelling of plasma reactions is im- 
plemented in the fluid context, using rate coefficients to 
determine the evolution of fluid plasma species. How- 
ever there is an implicit assumption in such fluid ap- 
proaches, namely that the underlying plasma distribu- 
tion is Maxwellian, with a wcU-deflned temperature; no 
macroscopic variables other than the density (the zeroth 
kinetic moment) enter such equations, and the tempera- 
ture therefore is an externally defined parameter of the 
problem. The rate coefficients are often presented in 
Arrhenius-type formulae that make the temperature de- 
pendence explicit, and the reaction set is then modelled 
by rate equations which are primarily driven by species 
number density variations within the Maxwellian con- 
text. For example 0, S [13, El > 



k ^ A{T/Tor exp[-Ea/{kBT)] 



(1) 



activation energy Ea, and reference temperature Tq; n is 
an additional parameter, theoretically predicted to take 
the value 1/2. An explicit example is that for the disso- 
ciative electron attachment to an ozone molecule, form- 
ing a negative atomic ion and a neutral molecule Ij, 1^ 



e -f O3 ^ O2 + O" 



k = 2.12 X 10-i5rg^i °6exp(-0.93/re)m3s 
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in which Tg is given in electronvolts. In this article we 
explore the possibility that the plasma reactions affect 
the plasma temperature (more correctly, the mean en- 
ergy per particle) in a feedback loop, simply because the 
distribution function of each species must be altered as 
a direct result of reactions converting one set of plasma 
species to another set. In order to describe this as a 
distribution function evolution, we require a gas-kinetic 
description based on the full reaction cross-sections, not 
to be confused with a chemical-kinetic one which uses 
the fluid-based rate constant, and which is an averaged 
treatment. In general, the rate constant k used in chem- 
ical kinetic descriptions is derived from the gas-kinetic 
reaction cross-section a by the following integration over 
velocity space v Q : 



k ^ 



J a{v)f{v, t)vdv 
Jf{v,t)dv 



(3) 



is the general Arrhenius form for a chemical process with 



where f{v,t) is the distribution function of the species 
(taken as Maxwellian) implicated in the reaction asso- 
ciated with the cross-section a, v is the velocity co- 
ordinate, and the integration is over the whole velocity 
space. The casting of the formal expression for k from 
Eq. ([3|) in the (modified) Arrhenius form of Eq. ([T]) is 
generally done by an empirical fit to experimental data, 
valid over a limited energy range. 

In this article we will concentrate on electron- 
moderated reactions, but the principle applies to all 



species. For instance, electron attachment and detach- 
ment can have a significant efi'ect on the electron distri- 
bution within a plasma; another example is recombina- 
tion, in which electrons and positive ions recombined to 
form neutrals, simultaneously depleting the positive and 
negative species whilst augmenting the neutrals. The 
fact that such reactions proceed according to a reaction 
cross-section in kinetic theory means that the evolution 
of the distribution function of each species can be ac- 
commodated in such a model; however, modelling in the 
fluid context with rate equations averages out distribu- 
tion function effects by requiring that the underlying dis- 
tribution is always Maxwellian for rate coefficients that 
are temperature dependent. The sections following this 
one describe the background kinetic theory, and show 
how the simplest possible model encompasses a temper- 
ature drift arising directly from the infiuence of the reac- 
tions on the form of the electron distribution function. 



BACKGROUND THEORY 

Consider a plasma in which binary interactions main- 
tain the equilibrium. Then the distribution function fg 
for a given species s obeys Boltzmann's equation. 



{v?){r, t) as moments of the distribution function: 



dt 



dr 



du 
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djs 
dt 



where ag is the acceleration, r is the spatial co-ordinate, 
and u is the velocity co-ordinate. The distribution func- 
tion is defined such that /(r, u, t)Au gives the number of 
particles in the velocity range u to u -\- du at the posi- 
tion r at time t. The collision term on the right gives the 
change in the distribution arising from binary collisions, 
and is given by Boltzmann's formula 1J| : 



dt 



X! j [fs{r,u',t)fj{r,u'^,t) - fs{r,u,t)fj{r,Uj,t)] 



\u — Uj \ ajgduj 



(5) 



where fj denotes the distribution function of the tar- 
get particles, and ' denotes post-collision quantities; ajs 
is the collision cross-section for the interaction between 
particles j and s. The term in square brackets shows 
how the distribution function is changed by the flux of 
particles scattered into the velocity range u to u + du 
minus those scattered out of that range, as a result of 
the interaction governed by the cross-section agj- Note 
that a kinetic plasma in equilibrium has zero collision 
term. The kinetic theory of plasmas will yield bulk (that 
is, macroscopic or fluid) quantities by integrating over 
all velocity space. For example, we can define the num- 
ber density n(r,t) and mean square speed per particle 



n(r,t) = J f{r,u,t)du 



{u^)(r,t) = / u^f{r,u,t)du 



J fir,u,t) du 



(6) 



(7) 



In a Maxwellian equilibrium, the mean energy per parti- 
cle defines the temperature T: 



1 / 2 

— miu 

2 ^ 



^NdksT 



(8) 



where m is the particle mass, and is the number of 
spatial degrees of freedom. Very often the electron distri- 
bution is not actually Maxwellian, but the concept of an 
effective temperature is deduced from the mean energy 
calculated from the second moment in this way. 

The Boltzmann collision integral Eq. ([5]) has restricted 
validity: it is only a good description of binary, uncor- 
related interactions which take place over length and 
time scales much shorter than any other intrinsic vari- 
ation in / itself (such as plasma inhomogeneity or ex- 
ternal forcing). Since charged particles in a fully ion- 
ized plasma tend to interact collectively via simultaneous 
long-range Coulomb forces arising from the disposition of 
other charges, the binary interaction is less appropriate, 
and the Boltzmann collision term is often replaced with 
the more sophisticated Fokker-Planck treatment. How- 
ever for a weakly ionized plasma in which such many- 
body Coulomb interactions do not dominate over binary 
encounters with neutral species or other charged parti- 
cles, the Boltzmann description is still very relevant. 



SIMPLE MODEL: WEAKLY IONIZED PLASMA 

Consider the simple case of a weakly ionized kinetic 
plasma, governed by the Boltzmann collision integral, in 
which a particular reaction suddenly starts in an unbal- 
anced way. For example, this could be an electronegative 
plasma in which electron attachment becomes the dom- 
inant driving term under a set of prevailing local condi- 
tions, or perhaps a plasma begins to quench as a result 
of recombination becoming important as a local energy 
density drops. In each scenario, the key issue is that the 
reaction is initiated in a plasma in such a way that it 
is predominantly one-sided, that is, it is far from equi- 
librium for that reaction. For the sake of being specific, 
consider a plasma in which the temperature has fallen 
sufficiently that recombination suddenly becomes a sig- 
nificant reaction. The sequence of events we imagine to 
proceed as follows: 

1. The plasma cools below some critical energy thresh- 
old such that the balance between recombination 
and ionization is significantly altered to favour the 
former 
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2. Binary recombination dominates, leading to a loss 
of electrons and positive ions according to the en- 
ergy dependence of the cross-section 

3. The distribution function evolves according to the 
Boltzmann equation, with the collision integral 
explicitly yielding corrections to the distribution 
function on timescales shorter thn the equilibration 
time 

4. The distribution function equilibrates to a new 
mean energy (or temperature, if Maxwellian) and 
the reactions proceed according to the new, evolved 
form of the distribution function that governs the 
remaining charged particles 

The significance of this description lies in the fact that 
the mean energy per particle (proportional to the tem- 
perature, if Maxwellian) is evolving at the same time as 
the reaction proceeds, since if we consider the correc- 
tions Sf to the distribution function at time t arising 
from the loss of particles, then the instantaneous distri- 
bution function will not take the same form as the initial 
one: the form of the collision integral Eq. ((5]) precludes 
this for all but pathological forms for a. Relaxation back 
to a Maxwellian (for example) for fewer particles may 
well be possible, but will necessarily change the temper- 
ature, since there has to be a redistribution of particles 
in velocity (and energy) space. 

In kinetic terms, this behaviour can be readily accom- 
modated (in principle); however, in the fluid context this 
is problematical, since there is an implicit assumption 
that the distribution function is always in its equilibrium 
form, that is, Maxwellian. Fluid models therefore use 
rate constants, rather than reaction cross-sections, and 
so the feedback loop that connects the progress of the 
reactions to the evolving form of the participating dis- 
tribution functions is broken. Non-equilibrium effects in 
plasma chemistry are known to be significant 15, 1^ and 
the reconstruction of effective reaction rates from gas- 
kinetic simulations can provide information inaccessible 
from only fluid models [17|. 

If as a direct result of kinetic effects either the temper- 
ature drifts, or if the distribution functions of the reac- 
tants become non-Maxwellian, then the assumed isother- 
mal rate coefficients may be inaccurate, leading to cumu- 
lative discrepancies in the relative numbers of species. 

The simplest case we can do for illustrative purposes is 
the 1-velocity dimension kinetic plasma in which a par- 
ticular reaction becomes predominantly one-way, leading 
to electron losses (that is, the electrons are recombin- 
ing with positive ions or perhaps are attaching to neu- 
trals). Let's assume that the electrons and ions (or neu- 
trals) (denoted by subscripts e, i and n respectively) are 
Maxwellian initially, with temperature Tq, though the 
analytical framework is more general than this. Since 
the reaction is essentially one-sided, meaning that elec- 



trons are much more likely to be captured (and lost to 
the distribution) than scattered into a different velocity 
range, we can approximate the collision term by retain- 
ing only the part that describes the loss of electrons from 
a velocity element: 



dfe 

dt 



fefs\vi 



ladv, 



(9) 



where s denotes the species with which the electron is 
interacting (s — i for recombination, s = n for electron 
attachment to neutral species n, etc), and we have sim- 
plified the notation by dropping functional arguments. 
The cross-section for the dominant process under consid- 
eration is denoted simply by a. Of course, should the 
reaction be creating more free electrons, then the sign of 
the right-hand side of Eq. ^ is reversed. Note that there 
is also a self-term for the electrons that will lead to a re- 
arrangement of the electron distribution without chang- 
ing the total number of electrons; the importance of this 
self-relaxation term has to be judged in terms of the rel- 
ative timescale for the competing reaction. We will first 
evaluate the effect of the electron loss on the electron dis- 
tribution function as a result of the continuing reaction, 
and then address how this may be modified by relaxation 
processes, given that electron-electron relaxation is best 
treated as a simultaneous, collective interaction rather 
than a sequence of binary interactions. 

In order to better illustrate the issue at the heart of 
this article, we will move from what has been gener- 
ally valid to a specific case of Maxwellian electrons. Of 
course, other driven equilibria are perfectly possible, and 
can be described by the same theoretical framework, but 
the integrals are known and tractable for the Maxwellian 
plasma, and so we will report the results in this case. 

Taking Maxwellians at a common temperature T ini- 
tially for the electrons and ions. 



as ^ ms/{2kBTo) 



Cs = nosVcWvr 



(10) 

(11) 
(12) 



in which s = i or e, and tiqs is the initial species number 
density. We are approximating the cross-section as a step 
function in energy, such that a = CTq for ^-U speeds less 
than some critical speed (to be defined later), and a = 
otherwise. (In practice, if the critical speed is greater 
than the thermal speed, there is a negligible contribu- 
tion to the integral of a non-zero cross-section at values 
higher than the critical speed for a Maxwellian distribu- 
tion; hence taking the cross-section as constant for all 
energies doesn't significantly change the results.) The 
simplicity of this appproximation is that all of the key 
integrals can be done analytically, illustrating the con- 
cept. The Appendix gives the full form of the integrals; 
to minimise the disruption to the flow of the argument 
we will simply state the analytical results here: 
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+aonQ^noeTT ^^/af^ja^ exp[-(Qi + ae)We] (14) 
dt 



(15) 



This yields the total rate of change of fe resulting from 
the recombination process, leading to the evolved distri- 
bution 



UAt) « /e(0) + At 



fdfeiO) 



V dt 

Cee"^"^' {1 - noaoAt [\ve\ erf(y^|i;e|) 



(16) 



In order to maintain a physically meaningful corrected 
distribution, we need /e(At) > 0, and so we must restrict 
the validity of pB]) to values of such that 



r^oo•oA^|^;e| < 1 



(17) 



bearing in mind that 3> ae, and so the value of the 
error function can be taken to be unity, and neglecting 
the last term on the right-hand side of Eq. (I16p . which 
is not significant at finite jwel- This restriction can be 
accommodated assuming that an energy (or indeed ve- 
locity) dependent cross-section will fall sharply with en- 
ergy (or speed) to become negligible in value outside a 
critical energy range, so that it closely resembles a step- 
function (for example, some electron attachment reac- 
tions are only significant below a critical threshold en- 
ergy); the formal statement of Eq. ([TT]) stems from the 
fact that we hold the cross-section to be constant for all 
speeds. If instead we rephrase this to be cr = (Tq for 
|ue| < l/(nocroAt); cr = otherwise, then we have a pre- 
liminary basis for defining the critical speed mentioned 
earlier, and therefore formalising the approach. 

We can now ask the question: is the evolved distribu- 
tion given in Eq. (jl6p a Maxwellian? The answer is no, 
since if we take the function 



0(a,j/) = e-ni-ay'/') 



(18) 



as a simplified, scaled approximation of /e(Ai) expressed 
in terms of energy, rather than speed, we can see that 



log0(a,y) ~ -y- a^/y 



(19) 



A true Maxwellian would have only the first term on 
the right-hand side in Eq. (fT9|l : there is an additional 
small, non-zero curvature correction to the straight line 
expected from the log-linear plot, characterising the drift 
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FIG. 1: Plot of the function (j){a,y) from (|18p as a function of 
normalised energy y for a = (pure Maxwellian), a = 0.1 and 
a = 0.2, illustrating the drift from the pure Maxwellian as a 
function of collision correction. Note that non-zero values of a 
still give the appearance of a straight-line solution in the log- 
linear plot, apart from significant curvature at low energies, 
though clearly the inferred temperature is lower than the true 
Maxwellian case. 



from a Maxwellian as a result of the unbalanced binary 
interaction. Figure [1] gives a simple illustration of this 
effect. Log-linear plots of the distribution function as a 
function of energy with different values of the collision 
correction show that the latter increases the apparent 
gradient of the plot, showing that the inferred 'temper- 
ature' (proportional to the reciprocal gradient) is lower 
than that associated with the pure Maxwellian. This 
change from a Maxwellian shouldn't be too surprising, 
given the nature of the binary collision term in Eq. ([9]): 
the distributions of the electrons and target particles are 
involved in a (relative) velocity-scaled convolution prod- 
uct, and only pathological cases could yield a gaussian 
form. 

In order to quantify the phenomenological conse- 
quences of this correction, we can calculate key macro- 
scopic observables. By integrating Eq. (fTH]) over electron 
velocity space we can evaluate the accompanying rate of 
change of electron number density: 



due 
~dt 



-1/2 



1/2 



(20) 



To get the change in the mean energy of the particles, we 
need to multiply Eq. p4p by Ug and integrate again over 
electron velocity space, resulting in the expression 



d{nevl) 1 _i/2 
at 2 



a. 



1/2 



2ai -\- Q!e 



(21) 
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Hence we can approximate the density and energy evo- 
lution in time Ai as follows: 



din 

ne{At) w noe+[-g^] At 



(22) 



noe - TT " ' noeUo^ao [ — — ) (23) 



1/2 



At (24) 



'2a 



X 



2^ 



ai + 



1/2 



Ai 



(25) 



so long as the time interval At is sufficiently short that 
the distribution function has not changed significantly, 
that is, the perturbation in particle number density must 
be small. 

In order to interpret the change in the mean energy as 
a change in temperature, we need to consider the equiv- 
alent MaxwcUian distribution that has the same particle 
number density and total energy as the collision-amended 
one. The temperature Ti of this equivalent Maxwellian 
can be given in terms of the original one Tq by consider- 
ing the ratio of the second moment to the zeroth of the 
post-collision distribution function, Taylor expanded in 
time using equations (p3|) and (|25|) . This yields 



Ti AT 

T^r = 1 - ^ = 1 - noi^roAt 
-to Jo 



OLi 



naeiai + ae) 



1/2 



(26) 



showing that the temperature drops as the interaction 
causes particles to be lost to the distribution, consistent 
with the simple picture presented in Figure [T] 

Now we can quantify the permitted time interval for 
this approximation to be valid: the second term in 
Eq. (I26|) must be small, yielding an appropriate value 
for At, the distribution function evolution time. If this 
At is similar to the equilibration time, then it is rea- 
sonable to assume that the plasma will relax to a new 
Maxwellian in Ai, giving a formal statistical basis for the 
evolving temperature. However, the relaxation might not 
be rapid in a weakly ionized plasma, since the electrons 
and ions are dominated by neutrals, and therefore the 
relatively rapid Coulomb interaction which would equi- 
librate a fully ionized plasma may in fact be dominated 
by electron-neutral collisions in such a strongly-coupled 
plasma (a plasma is weakly coupled if the kinetic energy 
of the particles greatly exceeds any local potential energy, 
and strongly coupled if the kinetic energy is dominated 
by local potentials). 

Hence there are several possibilities for the evolution 
of the plasma in time Ai: 



• the electrons and the positive ions could both relax 
to a new Maxwellian with a new temperature; 

• the electrons might relax to a Maxwellian, but the 
ions remain non-Maxwellian; 

• neither species relaxes to a Maxwellian 

Note that the relaxation of electrons must be via interac- 
tion with species other than the target one implicated in 
the reaction, since interaction with the latter is demon- 
strated to produce an evolution away from a Maxwellian. 
The most rapid relaxation process is self-relaxation, given 
the superior mobility of electrons compared to any other 
species present, and this process is the most likely to 
produce a Maxwellian form for the electrons. 

In the first two cases, Eq. ([26]) is an approximate 
Taylor-series expansion of the electron temperature evo- 
lution with time. In order to quantify the magnitude 
of this temperature change, it is necessary to specify a 
value for Ai. Since the assumption underlying Eq. (j26[) 
is that the electrons have relaxed to a Maxwellian, it is 
appropriate to take the time interval to be equal to the 
self-equilibration time r^" for electrons 3, Iq : 



{2kBT/mef^^ 

e^riQf, In A 
2TTelml 
erf (a:) — x erf '(a;) 
2^^ 



(27) 
(28) 
(29) 



where In A is the usual Coloumb logarithm. 

Recognising that Ue <^ ai, we can simplify the expres- 
sion for the fractional temperature change: 



AT 



To 4*(l)lnA 



ao 2^el{2kBTof 



3 X 10** X 



In A 



(30) 



where we have taken 44' (1) sa 1. The Coulomb loga- 
rithm is usually taken to lie in the range 2 < InA < 
20, with the upper limit associated with fully ionised, 
weakly coupled plasmas, and the lower one with partially 
ionised, strongly-coupled plasmas, though some caution 
has to be exercised here: the validity of such assump- 
tions has been challenged for plasmas that are not the 
idealised, Maxwellian ones for which the Fokker-Planck 
analysis holds perfectly (and therefore the concept of self- 
relaxation producing a In A term is consistent) [l^, |2^ , 
and in reality, the value of In A may be only half of the 
classical one. 

For a relevant specific example, let us return to the 
example of electron attachment to ozone, for which the 
rate constant in Arrhenius form was given in Eq. 
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The cross-section for this dissociative process has a peak 
at around 1 eV [l^ 21 1 (there is a resonance below 1 



eV which we shall disregard in the context of a plasma 
cooling to the process threshold at around 3 eV). The 
peak value of the cross-section (excluding the resonance) 
is ~ 3 X lO^^^m^. If we take ctq 2 x lO^^^m^, and take 
In A - 2, To - 3 X IQ-^K, then we find from Eq. ^ that 



AT 



3 X 10" 



(31) 



so that the fractional temperature change is under 0.1% 
on an electron equilibration time. A plasma with a signif- 
icant ozone fraction is clearly only partially ionized, and 
so the validity of using the equilibration time Eq. (f27|) is 
open to question. Nevertheless, this temperature drift is 
cumulative, since cooling the electrons maintains the like- 
lihood of the reaction. Thus in the fluid context, the min- 
imum temperature fluctuation must be around 1-2 orders 
of magnitude greater on the smallest fluid timescale, since 
the validity of the fluid approximation requires many col- 
lision times to elapse during the smallest resolvable fluid 
time. The significance of this lies in the evolution of the 
rate constant, which is a fluid quantity. The fractional 
rate of change of k in Arrhenius form Eq. ([T]) is 



(32) 



showing that the fractional drift in the rate constant k 
scales with the fractional temperature drift as a function 
of n and Ea, as well as T. Hence for a system of fluid 
species evolution equations, the drift in each rate con- 
stant can be different, leading to cumulative changes in 
the evolving species populations compared to those aris- 
ing under the usual assumption of uniform temperature. 

It is worth noting here that depending on the sign of 
the factor on the right-hand-side of Eq. ([5^ . the rate 
constant changes may be in or out of phase with the 
temperature drift. Hence it may be possible for the tem- 
perature to drift below a minimum threshold for a reac- 
tion, causing it to stop temporarily pending additional 
local heating. Such instabilities would only be apparent 
in the kinetic limit, but could have consequences for the 
overall efficacy of a plasma reaction process. 



ILLUSTRATING THE EFFECT OF KINETIC 
FEEDBACK ON THE FLUID MODEL 



In order to explore the consequences of Eq. f32l) , con- 
sider the simple chemical kinetics reaction 



A^ B 



(33) 



in which species A is modified by electron attachment or 
recombination to produce species B. Since the context 
here is the non-equilibrium evolution of the plasma, the 



reverse reaction has a negligible role to play. The rate 
equations for this reaction are 



?7<g — kTZ^TL 
Til — J\iT}jqT1ip^ 

— +knenp^ 



(34) 



where we are only taking equations for the one-sided re- 
action of Eq. p3p . with the electron number density, k 
the rate constant for the reaction and n^, the number 
densities of species A and B respectively. 

By eliminating Uj^ from the first two equations in 
Eq. ()34p . we arrive at a single equation for the evolution 
of the electron number density: 



(35) 



Integrating once, and neglecting the integration constant, 
we have 



which has solution 



kni = 



1 + / kit') At' 



(36) 



(37) 



where ne(0) = rieo and we have allowed the possibility of 
k being time-dependent, motivated by the explicit tem- 
perature dependence of rate equations shown in Eq. ([T]), 
and the fact that this temperature must drift in time, 
according to Eq. ([50)1 . 

Since the equation for is identical to that for rig, 
then species A evolves in time according to 



n^it) 



l + riAO / k{t')dt' 



(38) 



with the change in species B being determined via 

JT-A + = constant = ng, say (39) 

which is evident from the sum of the last two equations 
in Eq. Note that the evolution of in Eq. is 

constrained by Eq. 

Modelling the time dependence of k as 



fc = /cq + kit^ 



(40) 



where /cqj and p > are constants, we can contrast 
the electron number density at time t for the case ki — 
(rate constant has fixed value, implying isothermal evo- 
lution) with ki ^ (time-dependent rate constant for a 
non-isothermal evolution) by defining the ratio R{t) as 
follows: 



ne{t; ki ^ 0) 1 + neokot 



ne{t; ki = 0) 1 + Ueokot^ 

« l + /5(l-(f>); 

w 1/$; /3> 1 



(41) 

f3 < 1 (42) 
(43) 



6 



p=2 



FIG. 2: Plot of the function J?(/?) defined in Eq. (|4T| for 
p = 1 and p = 2, illustrating the increasing discrepancy in the 
electron number density with time for time-dependent rate 
coefficients compared to non-evolving rate constants. After 2 
characteristic times, the discrepancy is around 10% for p — 1, 
and around 20% for p = 2. 



where 



/3 = Ueokot 



is a characteristic time, and 



$ = 1 



ki f (3 



(44) 



(45) 



Hence depending on the elapsed time, the population 
of species can be significantly different if the plasma is 
not isothermal as the reaction progresses. An illustration 
of this effect is shown in [2l where the ratio R{/3) from 
Eq. (|4T|) is plotted for two cases: p = 1 and p = 2. In 
each case, the horizontal axis is the characteristic time /3, 
and (j) = 1 + O.la:^, simply for illustrative purposes. It is 
evident from the graphs that there can be a discrepancy 
in the electron number density of up to 20% after a couple 
of characteristic times. 

In Eq. ([55)1 we neglected an integration constant; 
restoring it yields a more general differential equation: 



with solutions in the form 



kCrip 



C 



1 - (1 - C/ueo) exp -C k{t')dt' 



(46) 



(47) 



in which the integration constant C is determined by ini- 
tial conditions. Although the form of Eq. (|47)) is more 
complex than Eq. ((37)) . the conclusions are nevertheless 
the same: the reaction yield as a function of time is differ- 
ent from the uniform case (although, of course, the final 
yield once the reactants have been exhausted is always 
the same). 



DISCUSSION 

We have shown here that for simple binary interactions 
leading to species abundance changes, gas kinetic theory 
shows that the distribution function of the interacting 
species is altered, on short tiniescales, by their participa- 
tion in the reactions. The result is that the mean energy 
associated with the reactants can change, influencing the 
progress of the reaction itself. For the idealised situa- 
tion in which the species start off as Maxwellians, this 
feedback mechanism changes the temperature of the re- 
actants (particularly the electrons) and can cause the re- 
actant yield to alter significantly. If such calculations are 
carried out in the fluid limit using rate equations and 
rate constants, the latter can be shown to be subject to 
measurable changes even in the fluid limit. 

In the zeroth-order fluid description, governed by rate 
equations, the species densities evolve without a self- 
consistent impact on the kinetic distribution of reactants 
(and therefore on the rate at which the reactions pro- 
ceed). This article illustrates the risk in such an ap- 
proach, but suggests a possible pragmatic way forward: 
rather than a full gas-kinetic simulation for at least the 
electrons, it may be possible to vary the rate coefficients 
in a way that reflects the underlying kinetic physics in 
the fluid limit. 

In a more realistic situation, there will be multiple re- 
actions underway simultaneously, each with a different 
set of cross-sections (in the kinetic limit) or rate con- 
stants (in the fluid limit) and each having a unique evo- 
lution under the changing distribution function. 

Reaction-rate uncertainties due to non-equilibrium ef- 
fects in plasmas have been postulated before [for example 
22| and demonstrated in experiments 23], and there are 
numerical simulations of reaction rates via kinetic PIC 
codes for RF plasmas that show the limitations of 
depending on the chemical kinetic fluid approximations. 

In practice, yield differences are apparent only for in- 
tegration times of a few /3; given that the typical elec- 
tron number density in discharges is 10^*~^^m~'^, and 
that rate constants tend to have values in the range 
10^^^ — 10~^^m^s~^, then signficant drifts will occur 
over periods longer than 100 ms. This is too long for 
electronegative- type instabilities Eq. [l|, 0, US] , but is less 
than the typical time for ozone destruction Eq. 25], for 
example. Hence ozone creation and destruction is a good 
test of the reasoning in this article, and there are several 
examples in the literature [1^ [l^, [S^, [111 which exper- 
iment and modelling are directly compared. The simula- 
tions are all based on the hydrodynamical rate equations, 
similar to Eq. I|34[) . though more extensive. In general, 
the qualitative agreement is very good, but the quanti- 
tative numerical predictions are not always reliable and 
accurate, with the calculated yields diverging from the 
experimental measurements significantly. Had the nu- 
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merical simulations all been exact, then there would be 
no room for speculation that drift in the rate expressions 
may be relevant. 

The simple fluid example given in this article is purely 
illustrative; in a real plasma there will be many such 
competing reactions. However, any temperature drift re- 
sulting from binary interactions will impact on species 
numbers across all such reactions. In the fluid approxi- 
mation, the consequences for the rate equations are sig- 
nificant due to the rate of change of any one species being 
proportional to the product of the participating species 
number densities and the effective rate constant at that 
time. Of course, if there are no relaxation processes avail- 
able to maintain the initial form of the distributions in 
all species, then even allowing for a time evolution of 
the rate constant will not be accurate enough to model 
reaction yields. 
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KEY INTEGRALS 
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The following integrals are central to the results of this 
article, and are gathered together here for convenience. 
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